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1. Introduction 

The strong cosmic censorship hypothesis says that the maximal Cauchy development 
of generic initial data for the Einstein equations should be inextendible. Since the maximal 
Cauchy development is the largest region of spacetime which is uniquely determined by 
initial data, this is essentially the statement that the time evolution of a spacetime can 
generically be fixed by giving initial data. The question of the validity of strong cosmic 
censorship is thus the question of predictability in general relativity. The above definition 
is in fact only one partial formalisation of the general idea of strong cosmic censorship 
and is based on the suggestion of Eardley and Moncrief[5] that the right way of looking at 
the problem is as a question about global existence of solutions of the Einstein equations. 
That this is so will be taken for granted in what follows. 

In any study of cosmic censorship where solutions of the Einstein equations are spec- 
ified by giving initial data the choice of matter model is very important. The opinion has 
been expressed in the past that it is best to restrict to vacuum spacetimes, at least in the 
beginning. If cosmological spacetimes are considered then the direct physical relevance of 
vacuum models (representing a completely empty universe) is questionable. Thus from 
the point of view of physics the best justification for studying vacuum cosmological models 
seems to be the idea that what is learned in the process may usefully be applied at some 
later date to non-vacuum spacetimes. In the situations considered in the following the 
vacuum solutions appear to play an exceptional role and so some caution is required when 
treating vacuum cosmological models as indicators of the behaviour to be expected in cases 
when matter is present. At this time there are unfortunately no results available on the 
global structure of cosmological spacetimes containing matter comparable in generality to 
those which have been obtained in the vacuum case[3,4,8,13]. 

As soon as non-vacuum spacetimes are considered, certain problems arise with cosmic 
censorship. Phenomena like shell-crossing and shocks lead to singularities (or at least to 
situations which in our present state of mathematical sophistication we are forced to con- 
sider as singularities in rigorous analytical investigations) which intuitively are due purely 
to the matter itself and have little to do with gravitation. These may obscure the gravita- 
tional singularities which are the true object of interest in the study of cosmic censorship. 
These difficulties have been discussed in detail in [20] where evidence is presented that a 
good candidate for a matter model which avoids these problems is the coUisionless gas, with 
the dynamics of matter being described by the Vlasov equation. In the asymptotically flat 
case a number of things have been proved about the properties of solutions of the Einstein 
equations coupled to matter of this type [14-17,19,21]. In the cosmological case much less is 
known although a general theorem concerning the long-time behaviour of a self-gravitating 
coUisionless gas has been proved in the context of Newtonian cosmology [18]. The aim of 
this paper is to begin investigating the global dynamics of cosmological solutions of the 
Vlasov-Einstein system. 

The simplest cosmological models are those which are spatially homogeneous and 
it is natural to look at these first. They fall into two broad classes, the Bianchi and 
Kantowski- Sachs models. The former can further be subdivided into Bianchi types I- 
IX. (For general information on this see [23].) There is a huge literature concerning the 
dynamics of homogeneous cosmological models which generally relies on one of the following 
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three methods: 

a) finding expUcit sohitions 

b) analysing a dynamical system in a small number of variables 

c) replacing the full equations by some approximate ones. 

The present work aims to obtain rigorous results concerning general solutions in a given 
symmetry class and so methods a) and c) are ruled out. Furthermore, even in the spatially 
homogeneous case the Vlasov-Einstein system involves infinitely many degrees of freedom 
(a function of the momentum variables) and so b) is also not an option here. A direct 
approach is necessary. The general idea is to first obtain the maximal Cauchy development 
by proving long-time existence in a constant mean curvature (CMC) slicing; in the homo- 
geneous case the CMC slices coincide with the hypersurfaces of homogeneity. Once this 
has been done inextendibility can be checked by determining whether incomplete causal 
geodesies run into a curvature singularity. 

The paper is organised as follows. Section 2 contains some generalities on Bianchi 
models and a proof that the occurrence of singularities in solutions of the Vlasov-Einstein 
system with a symmetry of this kind is controlled by the mean curvature of the hypersur- 
faces of homogeneity. A corresponding result for solutions with Kantowski- Sachs symmetry 
is obtained in section 3. These two sections also provide information on the qualitative 
behaviour of the mean curvature for the various symmetry classes. In section 4 it is shown 
that except in vacuum any singularity in a solution of the Vlasov-Einstein system with 
Kantowski- Sachs or Bianchi IX symmetry must be a curvature singularity, thus proving 
strong cosmic censorship in these two cases. It is also shown that the initial singularity 
is a curvature singularity for the other Bianchi types of Class A, once again under the 
restriction that the solution not be vacuum. To prove cosmic censorship for the remaining 
Bianchi types of Class A it would be sufficient to prove the completeness of future-directed 
causal geodesies. This is difficult and is carried out only for the simplest case, Bianchi 
type I, in the last section. 

2. Bianchi models 

Bianchi models are studied in this paper as a first step towards understanding more 
general cosmological solutions of the Vlasov-Einstein system and so it is appropriate to 
address the general question of which solutions of the Einstein equations should be regarded 
as cosmological. A conservative definition, and the one which will be adopted provisionally 
here, is that a cosmological spacetime is one which admits a compact partial Cauchy 
hypersurface. With this definition, the only Bianchi types which can occur for a spatially 
homogeneous cosmological model are I and IX. However this does not mean that the 
other Bianchi types are irrelevant in the present context because of the existence of locally 
spatially homogeneous cosmologies. 

Definition A spacetime is called locally spatially homogeneous if each point has a neigh- 
bourhood which is isometric to an open subset of a spatially homogeneous spacetime. 

With this definition, most Bianchi types can be realised in locally spatially homoge- 
neous spacetimes possessing a compact Cauchy hypersurface. Discussions of which Bianchi 
types these are can be found in [1] and [6]. It can be shown that the universal cover of 
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a locally spatially homogeneous spacetime is spatially homogeneous and since passing to 
the universal cover does not affect the dynamics it is enough to study the dynamics in the 
simply connected case. Then for Bianchi models the spatial manifold is in a natural way 
a simply connected 3-dimensional Lie group. The Kantowski-Sachs case, where passing to 
the universal cover brings no special advantage, is considered separately. 

Let G be a simply connected 3-dimensional Lie group and {e,} a left invariant frame 
on G. Denote the dual coframe by {e*}. Consider a spacetime with underlying manifold 
G X /, / an open interval, and a metric of the form —dt^ + gij{t)e^ . The initial 
value problem for the Vlasov- Einstein system will now be investigated in the case that the 
metric has this special form and the distribution function depends only on t and v% where 
V* denotes the spatial components of the momentum in the frame {cj}. Initial data will 
be given on the hypersurface G x {0}. The constraints are 

R - k'^hj + {g'^h.f = 167rToo, (2.1) 
V%j - Vjig^'^hm) = -SnToj. (2.2) 



The evolution equations are 



dtgij = -'2kij, (2.3) 
dthj = Rij + {g^'^kim)hj - 2kiik] - STTTij - ATTTooQij + 47T{g^'^Tim)9ij. (2.4) 

These equations are written using frame components. The second fundamental form is 
denoted by kij. R and Rij are the Ricci scalar and Ricci tensor of the three-dimensional 
metric respectively. The components of the energy- momentum tensor are denoted by Too, 
Toi and . The explicit form of the latter in terms of the phase space density / of particles 
is 

Too{t) = J f{t,v''){l+grsv''v'f/\AelgY'^dv^dv^dv\ (2.5) 
Toi{t) = J fit,v'')v^idetgy/^dv^dv^dv^ (2.6) 
(t) = J f{t,v^)viVj{l + grsv''v')-^/'^{dei gy/^dv^dv^dv^. (2.7) 



The Vlasov equation has the following form. 



df/dt + [2ky - (1 + grsv''v')-^/^-il^^v'^v'']df/dv' = 0. (2.8) 

The Ricci rotation coefficients 7*^ are defined in terms of the structure constants Cj^ of 
the Lie algebra of G by the following relations. 

^jk = \9^\^iik + Cjik + Ckij)-, (2.9) 
Cijk=9kiC\y (2.10) 
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To have a complete set of equations it is necessary to compute Rij in terms of gij but for 
the moment it is enough to know that Rij is of the form (det^)~"x (polynomial in gij 
andq,). 

A local existence theorem will now be proved under the assumption that initial data 
idij-i ^ij-i /°) given with f^{v) a function of compact support. The characteristics 
of (2.8) are the solutions V'^{s,t,v) of the equation 

dV'/ds = 2kiV^ - (1 + 9rsV'V')-^/''^i^^V'^V'' (2.11) 

with y*(t, t, f *) = f \ To control the determinant of the metric the following consequence 
of (2.3) will be used. 

d/dt{logdet g) = -2g'^kij. (2.12) 

The mean curvature g^-'kij will be of particular importance in the following and will be 
denoted from now on by H. Now define an iteration as follows. Let fo{t,v) — f^{v), 
goit) = g^ and kQ{t) = k^ . li fn, dn and k^ are given for some n determine V^+i by solving 
the characteristic equation (2.11) with g and k replaced by gn and Let 

fn+l{t,v) = f{Vn+l{Q,t,v)). (2.13) 

Next define an energy-momentum tensor Tn+i^ap by replacing / by fn+i and g by gn 
in (2.5)- (2.7). Now determine gn+i and kn+i by solving the linear ordinary differential 
equations which result when T^p, g and k arc replaced by T„+i^q,^, g^ and /c„ respectively 
in the right hand side of (2.3) and (2.4) and g and k by g-n+i and fcn+i on the left hand 
side. Let [0, T^+i) be the maximal interval on which gn+i is positive definite. By induction 
it can be shown that g^ and k^ are on their domains of definition. 

Let 1^1 be the maximum modulus of any component gij with a similar definition for 
kij . Now suppose that for all n < A?" — 1 the following bounds hold. 

\gn-g''\<Ai'i 
(dctgn)-^ <A2 \ (2.14) 
Ikn-k^'lKAs J 

Suppose further that \v\ < A^^ whenever fn{t,v) ^ 0. Here A\-A^ are positive constants 
which are for the moment arbitrary. The characteristic system (2.11) implies a bound of 
the form 

\v\ < Po + whenever /iv(t, v) i- 0. (2.15) 

where S4 depends only on A1-A4. As a consequence (2.5)-(2.7) imply a bound for TN,ai3 
depending only on A1-A4. Next (2.3) and (2.4) imply bounds of the form 



\kN - k^\ < Bst 



(2.16) 
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where Bi and S3 depend only on A1-A4. If A1-A4 (and hence Bi) are fixed then the 
inequahties (2.14) imply an inequality of the form 



(det^Tv)"' <S2 (2.17) 

whenever t <T and T is some positive time depending only on A1-A4. Now fix A1-A4 in 
such a way that A2 > {detg^)~^ and A4 > Pq. Next reduce the size of T if necessary so 
that BiT < Ai, B2< A2, B3T < A3 and Pq + B^T < A4. Then it follows that aU iterates 
exist on the interval [0, T) and that Qn and kn are bounded on that interval independently 
of n. 

Consider next the difference of successive iterates for n > 2. 



\9n+l{t) - gnit) \ + \kn+l{t) - kn{t)\ 



< 



C / \gn{s) - gn-l{s)\ + \kn{s) - A;n-l(s)| + ||/n+l(s) " fn{s)\\oods. 

Jo 



(2.18) 



For the difference of the characteristics there is the estimate 

\dVn+l/ds - dVn/ds\ < C{\Vn+l - K| + \gn " Qn-l] + \kn - /Cn-l|) (2.19) 

Let 

anit) = SUp{|K+l - Vn\{s,t,v) :0 < S <t,V E SUpp/n+l(t) USUpp/n(t)} gQ) 
+ \gn+l{t) - gn{t)\ + \kn+l{t) - kn{t)\ 

Then 

Wfn+lit) - /n(t)||oo < WfWc^anit). (2.21) 

Inequalities (2.19)-(2.21) imply that 

an{t) <C [ an{s) + an-i{s)ds (2.22) 
Jo 

Applying Gronwall's inequality to this gives 



Oin{t) < C / an-i{s)ds (2.23) 







It follows from this that an{t) < C"^~^||ci2||oo^"'~^/("' — 2)! so that {gn}, {kn} and {Vn} 
are Cauchy sequences on the time interval [0,T). Denote the limits of these sequences 
by g, k and V respectively. Going back to (2.3) and (2.4), we see that dtgn and dfkn are 
uniformly convergent. Thus /(t, v) = /'^(^(O, t, v)) is and ((/, /c, /) is a classical solution 
of equations (2.3)-(2.5) on the interval [0,T). If two solutions with the same initial data 
are given define a quantity a{t) in terms of their difference in the same way that a^it) was 
defined in terms of the difference of two iterates. Going through the same steps as above 
leads to an estimate of the form a{t) < C Jq a{s)ds. By Gronwall's inequality this means 
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that a{t) is zero and hence that the two solutions agree. Thus the solution which has been 
constructed is uniquely determined by the initial data. If the initial data satisfy (2.1) and 
(2.2) then so does the solution of (2.3)-(2.5). To show this define 

C = R- k'^kij + [g'^kijf - IGttToo (2.24) 
Ci = V'hj + SirToi (2.25) 

and note that (2.3)- (2.5) imply homogeneous first order ordinary differential equations for 
C and Ci. 

A close examination of the above proof shows that the size of T is only restricted by 
the quantities: 

\g%idetg')-\\k%\\f\\^ (2.26) 

and the diameter of supp/°. Thus if on some time interval [0,Ti) the quantities \g\, 
{detg)~^, ll/lloo and diam supp/ are bounded a solution exists on [t,t + e) for any 
t G [0,Ti) and some e independent of t. It can be concluded that the original solution can 
be extended to the larger interval [0, Ti + e). The following result has now been proved. 

Lemma 2.1 Let {g^, k^, f^) be an initial data set for equations (2.3), (2.4) and (2.8) which 
has Bianchi symmetry, satisfies the constraints (2.1) and (2.2) and is such that f^{v) is a 

function of compact support. Then there exists a unique corresponding maximal 
solution {g,k,f) on an interval [0,T*). If \g\, {detg)~^, \k\ and diam supp/ are bounded 
on [0,T*) then T* = oo. 

In order to investigate the question of global existence the mean curvature will be 
examined. 

dtH = R + H"^ - 127rToo + ^Trg'^Tij (2.27) 

This equation can usefully be combined with the Hamiltonian constraint (2.1). An equation 
which can be obtained in this way is 

dtH = k'^kij + 47r(Too + g^'Tij) (2.28) 

which shows that H is non-decreasing. Now suppose that H is bounded on [0, T*). Then 
from (2.28) 

H{t)>H{{))+ I {k'^kij){t')dt'. (2.29) 

^0 

Thus j^* ( k'^ kij ){t)dt< oo. 

At this point a short interlude on linear algebra is necessary. Let A and B he n x n 
symmetric matrices with A positive definite. For any matrix C the norm is defined as 
usual by 

||C||=sup{||Ca;||/H:x^O}. (2.30) 
It is also possible to define a relative norm by 

\\B\\a = sup{||Sa;||/||^a;|| ix^O}. (2.31) 
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It follows immediately that ||5|| < ||S||^||^||. The other fact which is needed here is that 



\B\\a < ^/tr{A-^BA-^B). 



(2.32) 



This can be proved by noting that there exists a basis {^i, ^2, ^3} with the property that 
Bzi = XiAzi for each i. The left hand side of (2.32) is then given by the maximum 
modulus of any A^. On the other hand, the eigenvalues of A~^BA~^B are {Af} so that 
the right hand side of (2.32) is equal to (Xli-^?)^^^- Let \\g\\ and ||A;|| denote the norms 
of the matrices with entries Qij and kij respectively. Let \\k\\g be the relative norm of the 
matrix with entries kij with respect to the matrix with entries gij. Then 



||^(t)||< 11^(0)11 +2 / \\k{s)\\ds 



< 11^(0)11 + 2 / ||/c(.)||,(,)||^(.)||d. 

^0 

< 11^(0)11 + 2 [\k^^h,y/\s)\\g{s)\\ds 

Jo 

Applying Gronwall's lemma to (2.33) gives 



(2.33) 



11^(011 < 11^(0)11 exp 



2 / {k'^kijy/^{s)ds 
. Jo 



(2.34) 



By what has been said above it can be concluded that if T* is finite then \g\ is bounded on 
[0, T*). Equation (2.12) shows that (detgf)"^ is bounded on that interval. It is known that 
if det g and its inverse are bounded the scalar curvature R is bounded from above [9]. Thus, 
from (2.1), k^^kij is bounded on [0,T*). Using the inequality \\k{s)\\ < (/c'-'/cij)^/^||fir(s)|| we 
see that |A;| is bounded. The boundedness of \g\ and {detg)~^ implies that g is uniformly 
positive definite on the interval of interest. Thus the solutions of the characteristic equation 
are bounded there. These statements add up to a stronger version of Lemma 2.1. 

Lemma 2.2 Suppose that the hypotheses of Lemma 2.1 hold. If H is bounded on [0,T*) 
then T* = 00. 

Lemma 2.2 shows that the long time behaviour of H gives useful information about 
the long time behavioiu' of the solution. Assume now without loss of generality that 
H < for t — 0. If this is not true for a given solution it may be arranged by doing 
the transformation 1 1— > —t. For all Bianchi types other than IX the scalar curvature R is 
non-positive [9]. Hence equation (2.1) implies that for those types H cannot become zero 
except in very special circumstances. More specifically, H = implies that Tqq = (so that 
the solution is vacuum), kij = and R = 0. Reference to [9] shows that for these Bianchi 
types the only cases where -R is possible are types I and VIIq. Moreover the metrics 
for which this condition is satisfied are flat. Thus the spacetime must be Minkowski space. 
To sum up: if the Bianchi type is not IX and the spacetime is not flat then iif < for all 
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time. Lemma 2.2 then implies that T* = oo. Next note that k'^^kij > l/3{g'^^kij)'^ so that 
(2.28) gives the inequahty 

dtH > \H^. (2.35) 

Comparing the solution of this equation with the solution of w = l/3tt^ with the same 
initial value shows that if H ever becomes positive then it diverges to -\-oo in finite time. 
Thus if it were known that H always became positive in type IX the ultimate fate of H 
would be known in all cases. That this is the case has been shown by Lin and Wald[10,ll] 
who proved under certain hypotheses (which are all satisfied by a solution of the Vlasov- 
Einstein system) that in type IX the mean curvature cannot remain negative for an infinite 
time. If it never becomes positive then, since it is non-decreasing, it must be identically 
zero for all t greater than some to- However it follows from the equations that in this case 
the 3-dimensional metric must be fiat, contradicting the fact that the Bianchi type is IX. 

S.Kantowski-Sachs models 

In terms of Gauss coordinates based on a hypersurface of homogeneity, the Kantowski- 
Sachs metrics are of the form: 

where the functions a and h depend only on t. They are spherically symmetric. Regular 
frames invariant under this symmetry do not exist and so an approach must be used which 
is different from that of the previous section. The Vlasov equation will be written in terms 
of the locally defined orthonormal frame: 

d/dt, a-^d/dr, b'^d/de, b'^ cscOd/dcj). (3.2) 

When written in terms of this frame the Vlasov equation for general functions in a geom- 
etry of the form (3.1) involves coefficients which are not regular functions on spacetime. 
However for distribution functions which are spherically symmetric these disappear. In 
the case of distribution functions / with full Kantowski- Sachs symmetry, i.e. spherically 
symmetric and independent of x, the equation takes the form: 

df/dt - a-^av^df/dv^ - b-^b{v^df /dv^ + v^df/dv^) = 0. (3.3) 

In an orthonormal frame the expressions for the components of the energy momentum 
tensor are standard. 

p = J /(l + |t;|2)l/2^^ldvW, 

p, = J f(v'fil + \v\Y^/^dv^dv^dv\ 

Here p is got by contracting T^p twice with eo and pi is got by contracting it twice 
with Cj. Equations (2.1)-(2.4) apply to the present situation provided the indices are 
now interpreted as coordinate indices. A local existence theorem for solutions of the 
Vlasov-Einstein system with Kantowski-Sachs symmetry can now be proved by following 



9 



the pattern of the proof of Lemma 2.1. The arguments leading up to Lemma 2.2 are also 

applicable in the present case. In fact the arguments are more transparent since the metric 
and second fundamental form are diagonal. The result of all this is the following: 

Lemma 3.1 Let {g^, /°) be an initial data set for equations (2.3), (2.4) and (2.8) which 
has Kantowski- Sachs symmetry, satisfies the constraints (2.1) and (2.2) and is such that 
f^{v) is a function of compact support. Then there exists a unique corresponding 
maximal solution {g,k,f) on an interval [0,T*). If H is bounded on [0,T*) then 
T* = oo. 

In order to determine the global behaviour of H it is enough, as in the Bianchi case, to 
find out whether in a solution where the function H is initially negative it will eventually 
be positive. An affirmative answer to this question has been obtained by Burnett [2]. As 
before the equation (2.35) then implies that H goes to infinity in finite time. Thus all 
Kantowski- Sachs solutions of the Vlasov- Einstein system exist only for a finite time in 
Gauss coordinates and the mean curvature of the hypersurfaces of homogeneity in such a 
solution takes on each real value precisely once. 



4. Curvature singularities 

For any (not necessarily spatially homogeneous) solution of the Vlasov- Einstein system 
the particle current density is defined by 

Ar« = - y fp"\ det ^^^g\^/^/podp^dp^dp^, (4.1) 

where p"^ denotes the coordinate components of the momentum. It satisfies 

Vc^AT" = 0. (4.2) 

Now suppose that some region of spacetime is foliated by compact spacelike hypersurfaces 
St with normal vector t"'. Then the quantity Jg^ N'^tadVg, where dVg is the volume 
element of the induced metric on the hypersurface, is time independent. In the case of 
a homogeneous spacetime it seems reasonable to hope that the integral can be dispensed 
with. For a Bianchi model takes the form: 

AtO = y"/(det^)VWd^;W. (4.3) 

Using this equation together with (2.8) and (2.12), the following identity can be derived: 

dt{{detgy/'N^) = -Cy^^Tokidetgy/^ (4.4) 

In the Class A models where, by definition, C*^ = the right hand side of (4.4) vanishes 
and (det g^^'^N^ is independent of time. For a Kantowski-Sachs spacetime (4.3) is replaced 

by 

N^ = J fdv^dv^dv^ (4.5) 
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and a short computation shows that once more {det g^^'^N^ is constant in time. 

Now consider a solution of the Vlasov-Einstein system of the type discussed in the 
previous two sections and suppose that for this solution H tends to infinity in finite time. 
Let T* be the time when H becomes infinite. The proof which was used to show that this 
happens in those cases was a comparison with an explicitly solvable ordinary differential 
equation. In fact more information can be extracted from this method. It can be seen 
that an inequality of the form H{t) > C{T* — t)~^ holds. This inequality implies that 
det g tends to zero as the singularity is approached. Hence for Kantowski-Sachs and Class 
A Bianchi models A^*^ tends to infinity. Comparing the definitions shows that T^^ > 
and so in these cases also diverges as the singularity is approached. Furthermore the 
scalar T^'^Tq,^ is greater than or equal to (T°°)^. The Einstein equations then imply that 
the curvature invariant C^^GafS diverges as the singularity is approached. By replacing t 
by —t this result can also be applied to what is, with the convention for the time direction 
chosen in section 2, the initial singularity. Thus the following result has been proved: 

Theorem 4.1 Let {g, k, /) be a non- vacuum locally spatially homogeneous solution of the 
Vlasov-Einstein system with Kantowski-Sachs or Bianchi Class A symmetry which is the 
maximal Cauchy development of data on a homogeneous hypersurface. Then every inex- 
tendible past-directed causal geodesic runs into a curvature singularity. If the symmetry 
type is Kantowski-Sachs or Bianchi IX then in addition every inextendible future-directed 
causal geodesic runs into a curvature singularity. 

This theorem shows that for Kantowski-Sachs and Bianchi IX models the maximal 
Cauchy development of any non-vacuum initial datum for the Vlasov-Einstein system is 
inextendible. Since non-vacuum initial data constitute a generic subset of all initial data 
for the Vlasov-Einstein system for any reasonable definition of genericity, this suffices 
to prove strong cosmic censorship within these two symmetry classes. It is, however, still 
interesting to ask what happens in the vacuum case. (When investigating cosmic censorship 
in inhomogeneous cosmological models it will be necessary to consider cases where there is 
matter somewhere and vacuum elsewhere.) The only vacuum Kantowski-Sachs model with 
a compact Cauchy hypersurface is obtained by identifying the part of the Schwarzschild 
solution inside the event horizon by means of the transformation t ^ t + C. Here t is the 
standard Schwarzschild time coordinate and C a positive constant. This spacetime has a 
curvature singularity in the future but in the past the Schwarzschild event horizon is now 
a Cauchy horizon, through which the spacetime can be extended to a region with closed 
timelike curves. Thus the vacuum solution behaves very differently from the solutions 
containing matter. In the Bianchi IX case it does not seem to be known in general which 
vacuum initial data have maximal Cauchy developments which are extendible. However 
there is an example where the spacetime is extendible in this way both in the past and 
the future, namely the Taub-NUT solution. Siklos[22] has shown that the Taub-NUT 
solution is the only only Bianchi type IX vacuum solution which can be extended to a 
homogeneous spacetime containing a Cauchy horizon and so it is tempting to conjecture 
that it is the only one which is extendible beyond the maximal Cauchy development of 
data on a homogeneous spacelike hypersurface. 
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5. Geodesic completeness 

To investigate the geodesic completeness of homogeneous spacetimes which expand 
indefinitely, consider first the following consequence of (2.27) and the Hamiltonian con- 
straint. 

dtH =1{R + H^ + Z0kij) + Airg'^Tij . (5.1) 

This section is concerned with spacetimes of Bianchi type I. In these spacetimes i? = 
and so (5.1) implies the inequality 

dtH > (5.2) 

thus improving on (2.35). Comparing the inequality (5.2) with the corresponding differ- 
ential equation gives: 

H{t)>-2{t + C)-^, (5.3) 
for some constant C. Hence it can be concluded from (2.12) that 

det^<Ct^ (5.4) 

It can also be concluded, using the Hamiltonian constraint, that 

Too < (l/47r)(t + C)-^ (5.5) 

A special feature of Bianchi I spacetimes is that the Vlasov equation can be solved ex- 
plicitly. Why this is so and why the same technique cannot be used for other Bianchi 
types is explained in [12]. For Bianchi type I the characteristic equation (2.11) implies 
that dVi/ds = 0. Thus if the distribution function is expressed in terms of Vi rather than 
the Vlasov equation just says that it is independent of time. Hence 

Too(t) = J f{v,){l + g'^%t)vrVs)'/\detg{t))-'/''dvidv2dv3. (5.6) 

Note that in this equation / does not depend on t. The matrix g'^^{t) is symmetric and so 
there must exist a basis which is orthonormal with repect to 5ij and consists of eigenvectors 
of this matrix. At any time there must be at least one of these which corresponds to 
the largest eigenvalue. Call it Wi. Then in the matrix notation introduced in section 2 
g'^^WiWj — \\g~^\\- Moreover 

g'^ViVj > g''^WrWs{5'hiWjf. (5.7) 
Assume that / is not identically zero and consider the expression 

FiVi) = J f{v,){d'H,yj)dv^dv2dv3. (5.8) 

If Ui is any vector which lies on the unit sphere in the sense that S^^yiyj = 1 then clearly 
F{yi) is a positive quantity which depends continuously on j/j. Hence it must have a 
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positive minimum m on the sphere and m depends only on /. Putting these facts together 
leads to the inequality 

Too{t)>m\\g-Y^\detg)-'/^ (5.9) 

Using (5.4) and (5.5) in this shows that H^"""^!! is bounded. Hence each component g'^^ is 
bounded and g^^ViVj is bounded along any causal geodesic. 

It is obvious that when a homogeneous spacetime described in Gauss coordinates 
exists for all t > then t must tend to infinity along any inextendible future-directed 
causal geodesic. Geodesic completeness can therefore be decided by looking at the relation 
between t and an affine parameter along such a geodesic. They are related by 

dr/dt = (/2 + g'^ViVj)-^/^ (5.10) 

where I is the length of the tangent vector to the geodesic, i.e. (— t'^ta)-^/^. Since g^^ViVj 
is bounded it follows that the right hand side of (5.10) can be bounded from below by a 
positive constant. Hence when (5.10) is integrated to get r the integral diverges as t — > oo 
and this is exactly what is needed to show future geodesic completeness. 

Theorem 5.1 Let {g, k, f) be a non- vacuum locally spatially homogeneous solution of the 
Vlasov-Einstein system with Bianchi I symmetry which is the maximal Cauchy develop- 
ment of data on a spacelike hypersurface. Then every inextendible future-directed causal 
geodesic is defined for arbitrarily large positive values of the affine parameter. 

Together with the results of section 4 this proves strong cosmic censorship in the 
class of Bianchi I solutions of the Vlasov-Einstein system. Once again it is interesting to 
ask what happens in the vacuum case. The vacuum equations for Bianchi I spacetimes 
can be solved explicitly to obtain the Kasner solutions [23]. It is straightforward to check 
that these spacetimes are future geodesically complete and that except in one special case 
(which is just a piece of Minkowski space in disguise) there is a curvature singularity in the 
past where the curvature invariant R'^^'^^ Ra/s-yS blows up. In the exceptional case there is a 
Cauchy horizon in the past. When the spatial topology is that of a torus the spacetime can 
be extended to a region containing closed timelike curves. The general vacuum solution 
has a very difi'erent behaviour at large times from the case where matter is present. It 
becomes more and more anisotropic as t — oo, contracting in one direction and expanding 
in two. The estimates for g'^^ obtained above show that when matter is present indefinite 
contraction in any one direction is impossible. In this sense the non-vacuum Bianchi type 
I solutions of the Vlasov-Einstein system resemble the Bianchi solutions with dust, which 
can be found explicitly [7]. 
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